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Abstract 

Top decays into a light Higgs boson and an up or charm quark can reach detectable 
levels in Standard Model extensions with two Higgs doublets or with new exotic 
quarks, and in the Minimal Supersymmetric Standard Model. Using both a standard 
and a neural network analysis we show that the CERN Large Hadron Collider 
will give 3<r evidence of decays with Br(i — » He) > 6.5 x 10~ 5 or set a limit 
Br(i — ► He) < 4.5 x 10~ 5 with a 95% confidence level if these decays are not 
observed. We also consider limits obtained from single top production associated 
with a neutral Higgs boson. 



PACS: 12. 15. Mm; 12.60.Fr; 14.65.Ha; 14.80.Cp 



1 Introduction 



The search for flavour-changing neutral (FCN) current processes both at high and low 
energies is one of the major tools to test the Standard Model (SM), with the potential 
for either discovering or putting stringent limits on new physics. In the SM, there are 
no gauge FCN couplings at tree level due to the GIM mechanism, and scalar couplings 
are also automatically flavour-diagonal provided only one Higgs doublet is introduced. 
However, there is no fundamental reason for having only one Higgs doublet, and for two 
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or more scalar doublets FCN couplings are generated at tree level unless an ad hoc discrete 
symmetry is imposed [1]. 

There are stringent limits on the FCN couplings between light quarks arising from the 
smallness of the K° — K°, B° — B° and D° — D° mass differences [2], as well as from 
the CP violation parameter e of the neutral kaon system. If there is no suppression of 
scalar FCN couplings, the masses of the neutral Higgs bosons have to be of the order of 1 
TeV. However, if FCN couplings between light quarks are suppressed, being for instance 
proportional to the quark masses [3,4] or to some combination of Cabibbo-Kobayashi- 
Maskawa (CKM) mixing angles [5,6], the lightest neutral Higgs may have a mass Mh of 
around 100 GeV. In this case, the effects of scalar FCN currents involving the top quark 
can be quite sizeable. 

The interactions among the top, a light quark q = u,c and a Higgs boson H will be 
described by the effective Lagrangian 

C = ^=g tq q(c v + c a j 5 )tH + h.c. (1) 



The axial and vector parts of the coupling are normalized to | c v \ 2 + | c a | 2 = 1 . The natural 
size of g tc in a two Higgs doublet model (2HDM) is g tc ~ 0.20, leading to large FCN 
current effects observable at the CERN Large Hadron Collider (LHC). The natural size 
of g tu ranges from 0.012 to 0.025, depending on the model considered. This is too small 
to be observed, but if this coupling is enhanced by a factor ~ 2, it could also be detected. 
Models with new heavy exotic quarks can also have large tree level Htq couplings, with 
the particularity that they are proportional to the FCN Ztq couplings [7]. Present limits 
on the latter [8] allow for g tu ~ 0.31 or g tc ~ 0.35, both observable at LHC. 

The large top quark mass opens the possibility of having relatively large effective vertices 
Htq induced at one loop by new physics. In the Minimal Supersymmetric Standard Model 
(MSSM) there are large regions of the parameter space where one can have g tc ~ 0.04 
[9,10] Q. In contrast, in the SM the effective Htc vertex is very small, g tc ~ 10~ 6 due 
to a strong GIM cancellation [11]. The large top quark mass also suggests that the top 
quark could play a fundamental role in the electroweak symmetry breaking mechanism 
[12,13] and in this case top FCN scalar couplings would be large, gt c ~ 0.3 — 0.4 [14] 
or 9tc ~ 0.5 — 2.9 [15], depending on the model considered. Hence the importance of 
measuring its couplings, in particular to the Higgs boson [16]. 

Present experimental bounds on top FCN scalar (and gauge) couplings are very weak. 
Some processes have been proposed to measure the Htq vertices, for instance tcv e v e and 
tce + e~ production at a linear collider with a center of mass (CM) energy of 1 — 2 TeV 
[17] and tc production at a muon collider [18] or at hadron colliders via gluon fusion 
gg — > H — > tc [14]. These processes can probe the FCN couplings of the heavier mass 
eigenstates, but do not provide useful limits for the light scalars. In this Letter we will 

1 The value of gt u is much smaller in principle. 
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show that top decays t — > if g at LHC provide the best limits on the couplings of a light 
scalar of a mass slightly above 100 GeV. The importance of this mass range stems from 
the fact that global fits to electroweak observables in the SM [19] seem to prefer M# ~ 100 
GeV, as well as the MSSM prediction that the mass of the lightest scalar has to be less 
than 130 GeV [20]. On the other hand, direct searches at the CERN e + e~ collider LEP 
at CM energies up to 202 GeV imply the bound Mh > 107.9 GeV for the mass of a SM 
Higgs, with a 95% confidence level (CL) [21]. For the lightest Higgs of the MSSM the 
bounds depend on the choice of parameters and typically are reduced to M H > 88.3 GeV. 

In the following we will perform a detailed analysis of the sensitivity of the LHC to FCN 
scalar couplings in top decays t — > Hq [22]. We will also study Ht production, which also 
serves to constrain the Htq vertex [23]. However, the bounds obtained from this process 
are much less restrictive. Finally, we will comment on the potential of the LHC to discover 
these signals at the rates predicted by some models. 



2 Limits on FCN couplings from top decays 

To obtain constraints on the coupling g tq we will first consider ti production, with the 
top quark decaying to W + b — ► Ivb, I = e, \i and the antitop decaying to Hq. With a 
good t identification the decays W + — > t + u can be included, improving the statistical 
significance of the signal. 

The decay mode of the Higgs boson depends strongly on its mass. For our evaluation, we 
take Mh = 120 GeV, well above the LEP limits, and decaying into bb pairs. In principle, 
there is no reason to assume that the Hbb coupling is suppressed. In a general 2HDM, 
both doublets couple to the b quark and, unless some fine-tuned cancellation occurs, the 
coupling of the physical light Higgs to the b quark is proportional to the b mass, which 
we take m& = 3.1 GeV at the scale \i = Mh- To be conservative, we assume that the 
Higgs couplings to the gauge bosons are not suppressed and are the same as in the SM. 
Thus for Mjj = 120 GeV H decays predominantly into bb pairs, with a branching ratio of 
0.7. For M H > 130 GeV, the W + W~ decay mode, with one of the W's off-shell, begins 
to dominate The signal is then lubbbj, with three b quarks in the final state. Clearly, 
b tagging is crucial to separate the signal from the backgrounds. We assume a b tagging 
efficiency of 50%, and a mistag rate of 1%, similar to the values recently achieved at SLD 
[24]. 

We evaluate the signal using the full gg, qq — > ti — > W + bHq — > Ivbbbq matrix elements 
with all intermediate particles off-shell. We calculate the matrix elements using HELAS 
[25] code generated by MadGraph [26] and modified to include the decay chain. For 

2 Of course, if the couplings to the gauge bosons are suppressed, for instance if H is a pseu- 
doscalar, the signal will get an enhancement factor of roughly 9/7 and the bb decay mode will 
dominate also for larger values of Mh ■ 
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definiteness we assume c v — 1, c a — 0, but we check that the difference when using other 
combinations is below 1%. We normalize the cross section with a K factor of 2.0 [27]. 



The main background to this signal comes from ti production with standard top and 
antitop decays, t — > W + b — > lub, i — > W~b — > jjfe. This process mimics the signal if one 
of the jets resulting from the W~ decay is mistagged as a b jet. The matrix elements for this 
process are calculated with HELAS and MadGraph in the same way as the signal, taking 
again K = 2.0. Another background to be considered is Wbbjj production, for which 
we use VECBOS [28] modified to include energy smearing and trigger and kinematical 
cuts. For this process we take K — 1.1 [29]. Wjjjj production is negligible after using 
b tagging. For the phase space integration we use MRST structure functions set A [30] 
with Q 2 = s. 

After generating signals and backgrounds, we simulate the detector resolution effects with 
a Gaussian smearing of the energies of jets j and charged leptons I [31], 

AEl 10 % noo/ AE3 50 % oo/ 

© 0.3% , — — = -= © 3% , (2) 



E l JW ' ' Ei 



where the energies are in GeV and the two terms are added in quadrature. We then apply 
detector cuts on transverse momenta pr, pseudorapidities i] and distances in (77, <p) space 
AR. 

p l T > 15 GeV , j? T > 20 GeV 
IV I < 2.5, AR jjtlj >0A. (3) 



For the events to be triggered, we require both the signal and background to fulfill at 
least one of the trigger conditions [32]. These conditions are different in the low 10 fb _1 
(L) and high 100 fb _1 (H) luminosity Runs. For our processes, the relevant conditions are 

• one jet with p T > 180 GeV (L) / 290 GeV (H), 

• three jets with p T > 75 GeV (L) / 130 GeV (H), 

• one charged lepton with p T > 20 GeV (L) / 30 GeV (H), 

• missing energy E T > 50 GeV (L) / 100 GeV (H) and one jet with p T > 50 GeV (L) / 
100 GeV (H). 

We reconstruct the signal as follows. There are three tagged b jets bis i n the final state 
and a non-6 jet j. There are three possible pairs 6j, bj, with invariant masses M^., one of 
which results from the decay of the Higgs boson and has an invariant mass close to M#. 
The invariant mass of this pair and the jet j, M b . bjj , is also close to the top mass. Hence 
we choose the pair b iy bj which minimizes (M bib . — M H ) + (M bib .j — m t ), defining the Higgs 
reconstructed mass as M|f c = M bibj and the antitop reconstructed mass as m r f c = M bibj j. 
In Figs 1, 2 we plot the kinematical distributions of these variables for the signal and 
for the ti background. The mf c distribution is slightly broader for ti due to the signal 
reconstruction method. 
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Fig. 1. Reconstructed Higgs mass Mjf c distribution before kinematical cuts for the signal and 
ti background in LHC Run L. We use gt q = 0.2. 
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Fig. 2. Reconstructed antitop mass mf c distribution before kinematical cuts for the signal and 
ti background in LHC Run L. We use gt q = 0.2. 

The remaining jet bk and the charged lepton result from the decay of the top quark. 
To reconstruct its mass, we make the hypothesis that all missing energy comes from 
a single neutrino with p u = (E u ,p/p ,p v L ), and p/p the missing transverse momentum. 
Using (p l + p u ) 2 = Myy we find two solutions for p u ', and we choose that one making 
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Fig. 3. Reconstructed top mass m\ cc distribution before kinematical cuts for the signal and ti 
background in LHC Run L. We use gt q = 0.2. 

the reconstructed top mass m\ cc = \J {p l + p u + p bk ) 2 closest to m t . In Fig. 3 we plot the 
kinematical distribution of this variable for the signal and ti background. 

The reconstruction of the signal and the requirement mj ec ,mp c ~ m t , M v ^ c ~ Mh are 
sufficient to eliminate the Wbbjj background, but hardly affect ti, as can be seen in 
Figs. 1-3. To improve the signal to background ratio we reject events when one pair of 
jets, bij or bjj, seems to be the product of the hadronic decay of a W boson, as happens 
for the ti background. We define the reconstructed W mass as the invariant mass M b .j, 
M bjj closest to M w (see Fig. 4) and impose a veto cut on events with M^f ~ M w . 

The complete set of kinematical cuts for both Runs is summarized in Table 1. We also 
require a large transverse energy H T . 



Variable 


Cut 


Mf c 


110-130 


m\ ec 


160-185 


m\ cc 


160-190 


1V1 W 


< 65 or > 110 


Hx 


> 240 



Table 1 

Standard kinematical cuts for the Ivbbbj signal. The masses and the energy are in GeV. 
Alternatively, we can use an artificial neural network (ANN) as a classifier to distinguish 
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Fig. 4. Reconstructed W mass M^ c distribution before kinematical cuts for the signal and tt 
background in LHC Run L. We use gt q = 0.2. 

between the signal and the tt background. In addition to Mj| c , mf c , m r t cc and M^> c , we 
consider pff, the transverse momentum of the reconstructed Higgs boson, p™ ax and p™ m , 
the maximum and minimum transverse momentum of the jets, and p^ max , the maximum 
transverse momentum of the two fo's which reconstruct the Higgs boson. The kinematical 
distributions of these variables are shown in Figs. 5-8. 

We observe that these variables are clearly not suitable to perform kinematical cuts on 
them. However, their inclusion as inputs to an ANN greatly improves its ability to classify 
signal and background events. We have not found any improvement adding other variables, 
and in some cases the results are worse. 

To construct the ANN we use JETNET 3.5 [33]. We find convenient using a three-layer 
topology of 8 input nodes, 12 hidden nodes and 1 output node. However, we do not 
claim that either the choice of variables or the network topology are the best ones. The 
performance of the ANN is not very sensitive to the number of hidden nodes. The input 
variables are normalized to lie approximately in the interval [—1,1] to reduce the training 
time. The network output r is set to one for the signal and zero for tt. We train the 
network using two sets with roughly 40000 signal and 40000 background events and the 
standard backpropagation algorithm. To test the ability of the network to classify never 
seen data we use two other sets with the same size. 

One frequent problem is overtraining. To avoid it, we keep training while the network 
error evaluated on the test sets, e[ \ decreases for the successive training epochs N. 
When it begins to increase, we stop training. To avoid fluctuations, we keep track of E[ N ^ 
for previous epochs. When we reach an epoch N Q when £]^ No ~ 100 ^ j s smaller than all the 
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Fig. 5. Higgs transverse momentum pj. distribution before kinematical cuts for the signal and 
tt background in LHC Run L. We use gt q = 0.2. 
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Fig. 6. p™ ax distribution before kinematical cuts for the signal and tt background in LHC Run 
L. We use g tq = 0.2. 

following values, we stop training and start again until epoch N — 100. 

We repeat the same procedure using four other different training and test sets of the same 
size, and find that the results are fairly stable. The final training is done including the 4 
training and 4 test sets. At the end of the training, most of the signal and background 
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Fig. 7. p™ m distribution before kinematical cuts for the signal and tt background in LHC Run 
L. We use gt q = 0.2. 
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Fig. 8. p^ max distribution before kinematical cuts for the signal and tt background in LHC Run 
L. We use g tq = 0.2. 

events in the test samples concentrate in very narrow intervals [0.99,1] and [0,0.01], 
respectively. 



To classify signal and background we use two other signal and background sets with never 
seen before data, 50 and 100 times larger, respectively, than the ones used in training, 
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to avoid possible statistical fluctuations. To separate signal from background we simply 
require r > 0.998. This maintains 29% of the signal while it rejects 99.9% of the tt back- 
ground. It is not necessary to train another ANN to distinguish the Wbbjj background, 
the same network with r > 0.998 completely eliminates it. In Table 2 we collect the 
number of events without cuts, with the standard cuts in Table 1 and with the ANN cut 
r > 0.998, for LHC Runs L and H. We normalize the signal to g tq = 0.2. 







Run L 






Run H 






before 
cuts 


standard 
cuts 


ANN 
cuts 


before 
cut 


standard 
cuts 


ANN 
cut 


Signal 


267 


98.2 


76.2 


2150 


797 


614 


tt 


7186 


33.2 


10.0 


58230 


270 


80 


Wbbjj 


77 


0.3 


0.1 


644 


2.2 


1.0 



Table 2 

Number of Ivbbbj events before and after kinematical cuts for the signal and backgrounds. We 
use gtq = 0.2. 

Although for definiteness we train the ANN with a vector FCN coupling (c v — 1, c a — 0), 
we check that the same ANN correctly recognizes the signal events when we consider an 
axial, left- or right-handed FCN coupling. The differences between the cross sections are 
in all cases below 1%, and comparable to the statistical Monte Carlo uncertainty. The 
differences between the cross sections after standard kinematical cuts are also smaller 
than 1%. 

To derive upper bounds on the FCN couplings we follow the Feldman- Cousins construction 
[34]. For the numerical evaluation of the confidence intervals we use the PCI package [35]. 
If there is no evidence of this process, i. e., if the number of events observed n is equal 
to the expected background n b , we obtain independent limits on the couplings g tu , g tc . 
Using the standard cuts we obtain g tq < 0.071 at Run L and g tq < 0.041 at Run H. Using 
the ANN cut improves these limits, gt q < 0.064 at Run L and g tq < 0.035 at Run H. 
All bounds are calculated with a 95% CL. Although the improvement may not seem very 
significant, in Section 4 we will see that the luminosity required to see a positive signal is 
reduced to one half with the help of the neural network. 

It is worth to give here also the results for the Fermilab Tevatron. This process at Run 
II with a CM energy of 2 TeV and a luminosity of 20 fb _1 gives only g tq < 0.38 (using 
standard kinematical cuts) due to the smaller statistics available. 



3 Limits on FCN couplings from Ht production 

The existence of top FCN Higgs couplings leads to Ht production through the diagrams 
in Fig. 9. This process has a negligible cross section in the SM [36]. As in the previous 
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section, we consider only the leptonic decays of the top quark, and the signal is lubbb. 
Note that Ht production is completely analogous to Zt production via Ztq anomalous 
couplings with Z — > bb [37], and has the same large backgrounds as this decay mode. The 
most dangerous is tt production, with t — > W + b — > lub, t — > W~b — > jjb, when one of the 
jets resulting from the W~ decay is missed by the detector and the other is mistagged as 
a b jet. Another small background is Wbbj production. 




Fig. 9. Feynman diagrams for gq — ► Ht via FCN couplings. 



The Ht process has a smaller cross section than tt with subsequent decay t — > Hq for the 
same values of gt q , but also receives an additional contribution from the latter when q is 
missed by the detector. This correction is not very significant for the analogous case of 
Zt production, about +10% for the up quark and +50% for the charm, but in this case 
with Mu > Mz the additional contribution becomes more important. In our calculations 
we will then take into account both processes. The signals and backgrounds are generated 
as for the Ivbbbj case. To estimate the number of events in which a jet is missed by the 
detector, we consider that a jet is missed when it has a low transverse momentum jP T < 20 
GeV or a very large pseudorapidity if > 3. 

The signal is reconstructed choosing first the pair of b jets whose invariant mass is closest 
to Mjj. The remaining b jet is assigned to the decay of the top quark, whose mass is 
reconstructed as in the previous Section. To improve the signal to background ratio we 
use the kinematical cuts shown in Table 3. 



Variable 


Cut 




115-125 


m rec 


160-190 




> 220 



Table 3 

Kinematical cuts for the lubbb signal. The masses and the energy are in GeV. 

In this case we cannot reduce the background imposing a veto cut to reject events in 
which a pair of jets has an invariant mass around Mw, because one of these jets is missed 
by the detector. Using an ANN does not improve the situation. The number of events 
for the gu — > Ht and gc — > Ht signals, the additional contribution tt — > Ht(q) and their 
backgrounds is collected in Table 4. 

The limits obtained from the Ht process are much less restrictive than from top decays. 
If no signal is observed, we obtain from Table 4 g tu < 0.13, g tc < 0.19 at Run L and 
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Run L Run H 





before 
cuts 


after 
cuts 


before 
cuts 


after 
cuts 


gu — ► -fft 


57.2 


35.1 


470 


293 


c/c -► Ht 


11.9 


6.8 


95 


56 


tt -► fTt(g) 


42.1 


22.0 


170 


186 


tt 


1373 


143 


11260 


1190 


Wbbj 


77.2 


1.6 


619 


12.7 



Table 4 

Number of Ivbbb events before and after the kinematical cuts in Table 3 for the Ht signal, the 
additional contribution tt — > Ht{q) and their backgrounds tt and Wbbj. We use gtq = 0.2. 

9tu < 0.076, g tc < 0.11 at Run H. The same is true for Tevatron Run II, where the bounds 
obtained are very poor, g tu < 0.40, g tc < 1.2 



4 LHC discovery potential 

In the previous Sections we have considered that the signal is not seen, and have obtained 
95% CL upper bounds on the size of the FCN couplings. These bounds imply the limits 
in the branching ratios of t — > Hq in Table 5. 

LHC Run L LHC Run H Tevatron Run II 

Signal Bi(t->Hu) Br(t -> He) Br (t -> Hu) Br (t -> He) Br(t ->■ Hu) Br (t -> He) 

Ivbbbj 1.4 x 10~ 4 1.4 x 10" 4 4.5 x 10~ 5 4.5 x 10" 5 5.2 x 10~ 3 5.2 x 10" 3 

Ivbbb 6.6 x 10- 4 1.3 x lO" 3 2.2 x 10" 4 4.3 x 10~ 4 5.9 x 10~ 3 5.3 x 10" 2 
Table 5 

95% CL upper limits on branching fractions of anomalous top decays obtained from the two 
processes analyzed. We also include the results for Tevatron Run II for comparison. 

The limits obtained for different values of M# are very similar. For a lighter Higgs, the 
larger phase space (in the case of top decays) or the larger structure functions (in Ht 
production) increase the signal cross section and hence allow a better determination of 
g tq . This improvement disappears when expressing the limit in terms of a branching ratio. 
After repeating the analysis using M# = 110, Br(H — > bb) = 0.8, we find limits 3% 
smaller from top decays and 5% smaller from Ht production. A heavier Higgs has the 
opposite behaviour, with the disadvantage of a smaller branching ratio H — > bb. With 
M H = 130, Bv(H — * bb) = 0.45 the limits from top decays and Ht production are 1.5 and 
1.3 times larger, respectively. 

Let us turn now to the possibility of observing these anomalous decays. The most common 
criterion to estimate the discovery potential is to use the discovery significance defined 
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as the ratio of the expected signal divided by the expected background fluctuation, s = 
n-s/ \/ni>- This is adequate when n?, > 5 and the Poisson distribution of the background 
can be approximated by a Gaussian with standard deviation ^fn^. One has "evidence" 
(3<t) when s > 3, and "discovery" (5 a) occurs with s > 5. 

In a general 2HDM, the most common ansatz for the FCN couplings is to assume that they 
scale with the quark masses. Under this assumption, the expected size of these couplings 
is g tc ~ ^m t m c / M w = 0.20, g tu ~ ^m t m u / M w = 0.012. An Htc coupling of this size 
(leading to Br(t — > He) = 1.5 x 10~ 3 ) will easily be detected and precisely measured. In top 
decays, at Run H we have s = 48 using standard cuts and s = 68 using the network cut. 
The benefit of the ANN is clear: the luminosity required to discover a signal is reduced to 
one half. In Ht production, the significance is much smaller, s — 7. The minimum coupling 
which can be discovered with 5 a is g tc = 0.054 (Br(t — > He) = 1.1 x 10~ 4 ), and the 
minimum coupling which can be seen with 3 a is gt c = 0.042 (Br(t — > He) = 6.5 x 10 -5 ). 
An Htu coupling g tu = 0.012 is too small to be seen, because in top decays we have only 
s = 0.23, and in Ht production s = 0.046. Note that for a 2HDM this ansatz implies 
that there are other interesting FCN current processes, for instance H — > /it at a muon 
collider [38]. 

Another possibility is to consider that the FCN couplings between two quarks are propor- 
tional to some CKM mixing angles involving these quarks. Rephasing invariance implies 
that these couplings are of the form g qq i = WqiV*,^ with i — d,s,b [5,6]. Choosing i = b, 
the hierarchy among the mixing angles guarantees a small value for g uc = XV u bV* b . With 
A = 5 we recover the values of g uc and g tc of the previous ansatz, but g tu is a factor of two 
larger, g tu = 0.025. An Htu coupling of this size is now near the detectable level: s — 1.1 
in top decays and in s = 0.22 in Ht production, and if it is enhanced by a factor ~ 2 it 
would also be detected. 

In models with exotic quarks the top FCN scalar couplings can be large, g tu ~ 0.31 or 
g tc ~ 0.35 (but not both simultaneously). These are the only models where the top can 
naturally have a large mixing with the up quark. However, the effects of top mixing in 
these models are most clearly seen with the appearance of tree level Ztq couplings, leading 
to top anomalous decays t — > Zq [39] and Zt production [37]. If these processes are not 
observed at LHC, the Htq couplings (proportional to the FCN Ztq couplings) must be 
small, g tu < 0.024, g tc < 0.030, below the detectable level. 

To analyze the case of top decays t — > Hq in the MSSM it is necessary to recalculate 
the signal and backgrounds and train again the ANN with the new data. The reason is 
that in the region of parameters of interest the Higgs width is much larger due to the 
enhanced Hbb coupling, and the M^ c distribution is broader in principle. For consistency 
we will use the parameters of Ref. [10]: tan/? = 35, a pseudoscalar mass M A = 120 
GeV and c v = c a = l/\/2. With these parameters, the Higgs mass is in the interval 
Mu — 105 — 120 GeV [40]. The kinematics is fairly independent of the particular value 
of Mu used, as we have shown above, the only difference being the Higgs decay rate to 
bb which varies slightly in this mass range. For our analysis we take Mu = 110 GeV, 
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Br(H — > bb) = 0.9 [41]. Our limits on Br(£ — > Hq) are proportional to 1/Br(i7 — > bb). To 
separate signal from background we select events with network output r > 0.994. 

There are large regions in the MSSM parameter space where Br(£ — > He) is greater than 
1.1 x 10~ 4 and even reaches 4 x 10~ 4 [10]. These decays can be discovered at LHC Run H 
with 5cr — 18 a. Also, a rate Br(t — > He) = 6.7 x 10 -5 is observable with 3 a. The decays 
t — ► Hu have a too small rate to be seen. 

On the other hand, if no signal is observed we obtain a limit g tq < 0.054 at Run L and 
9tq < 0.031 at Run H. Note that although these limits seem better than those obtained 
for Mh = 120, when translated into branching ratios they give very similar bounds, 
Br(t -> Hq) < 1.4 x 10~ 4 at Run L and Br(t -> ffg) < 4.7 x 10~ 5 at Run H. 
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